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Abstract. We give operational formulae of Burchnall type involving complex Hermite 
polynomials. Short proofs of some known formulae are given and new results involving 
these polynomials, including Nielsen's identities and Runge addition formula, are derived. 
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1. Introduction 

Operational formulae, involving special orthogonal polynomials, appear as useful tools 
to obtain some new results or to give short proofs of known formulae, see Burchnall [HE], 
Carlitz [6], Al-Salam [2J, Gould-Hopper [12J, Chatterjea [7], Singh [21], Karande and 
Thakare [18], Al-Salam and Ismail [3], Patil and Thakare [2T], Chatterjea and Srivastava 
[H], Purohit and Raina [22], Kaanoglu [16] among others. As basic example, one cites the 
classical real Hermite polynomials ([T3J, [201 I2S1 127]): 

H m {x) = (-l) m e* 2 D m (e-* 2 ); D := A 

which are extensively studied and have found wide application in various branches of 
mathematics, physics and technology. They possess a complete and rich list of remarkably 
interesting properties and some of them can be obtained using the Burchnall's operational 
formula (0]), to wit 

Indeed, it can be employed to give a direct and simple proof of the Runge addition formula 

m, my- 

as well as of the quadratic recurrence formula (Nielsen's identity [TU]): 

min(m,n) , R , , „ , , 

H m+n {x)-m\n\ 2. k \ {m-n)\{n-k)V ^ 



An extension of (1.1) is given by Gould and Hopper in [T2] for the generalized Hermite 
polynomials 

Hl(x,a,p) := {-l) m x- a e pxl D m (x a e~ vx " 
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As an interesting extension of the real Hermite polynomials H m (x) are the complex 
Hermite polynomials H p>q (z,z), for z = x + iy G C; x, y e R. They were considered by 
Ito in [J3] in the context of complex Markov process, 

ff M (v) = (-lfV^(e-j, 

where d z and c\ stand for 

d 1 ( d d \ , 9 l/<9 9 



^ z dz 2 y<9x <9x J ^ <3z dz 2 \dx ~^^dxj ' 

For the unity of the formulation, we define trivially H m (z,z) = whenever p < or 
g < 0. This class of polynomials is connected to H m {x) by 

and constitutes an orthogonal basis of the Hilbert space L 2 (C; e~ zz dxdy) (see [HI ITT]): 

/ H p ^(z,z)H m ^ n (z,z)e~ zz dxdy = 8 prn 5 qn p\q\'K. 
Jc 

Moreover, they appear naturally when investigating spectral properties of some second 
order differential operators of Laplacian type ([221 HOI ESI ED]) and can be connected to 
some classes of generalized Bargmann spaces. Several interesting features of H Ptq (z,z) 
in connection of coherent states theory have been studied recently (see for examples 

[DEI EH]). 

In course of our investigation, we establish new operational formulae involving H Ptq (z, z) 
analogous to the Burchnall one (1.1). In addition, we will use them in a simple way 
to obtain new properties satisfied by these polynomials. Mainely, we are interested in 
Nielsen's identities, generating functions and Runge addition formula. 

2. Burchnall's operational formula for H p>q (z,z) 

In order to obtain some operational formulae for H p q (z,z) of Burchnall type, we begin 
by noting first that the complex Hermite polynomials 



H m (z,z) = {-\r^m (e- z ~ z ) (2.1) 

= {-\) v e zz dl (z q e~ zz ) (2.2) 

= {-l) q e zz d q z (z v e~ zz ) , (2.3) 
can be rewritten in the following equivalent form 

H m (z,z) = (-d^ + z) p (z q ) (2.4) 

= (-d z + z) q (zn, (2.5) 
which appears as special cases of the following 

Lemma 2.1. For every sufficiently differentiable function f , we have 

(-flSr + zf (f) = (-l)V^ (e--/) (2.6) 

(~d z + z) q {f) = (-l)V^f (e-*V) • (2-7) 



Proof. We note first that (2.6) can be handled in a similar way as down below for (2.7). 
Indeed, direct computation yields (— d z + z) (e zz g^j = e zz (—d z g) . Next, by induction, we 
get 

{-d z + zf (e zz g) = {-iye z ~ z dl{g) 



which infers (2.7) for g = e zz f. 
Consequently, it follows 



□ 



Proposition 2.2. For given positive integers p, q and every sufficiently differentiate func- 
tion f , we have 



k=0 



k\ (q-k)\ 



(2.8) 
(2.9) 



Proof. Application of Leibnitz formula to the product function e zz f in (2.7) gives 

g 



'-d z + zY{f) = {-Ife^Y. \lW k (e-«) d k (f) 
A (-D fe z q - k 



k\ Uj-k)l 



d k z (f). 



Now by writing / as / = z p (z p f) in (2.7) and next applying Leibnitz formula, we get 
(~d z + zf (/) = (-1) V* £ ( q \ dl k (fe- a ) d k (z~ p f) 



k=0 



□ 



As immediate consequence of (2.8) when taking / = z p , we get the explicit expansion 



of the complex Hermite polynomial, 



min(p,q) f_\\k -^f-fe z p~k 



i/ p „( 2 ,r)= P !„! g H {q _ ky(p _ k y 



(2.10) 



The expression ( |2.9 ) for the complex Hermite polynomial can be considered as an 
analogue of (1.1). However, below we obtain a more appropriate Burchnall's operational 
formula involving H p ^ q (z,z). 

Proposition 2.3. For given positive integers p, q and every sufficiently differentiable func- 
tion f , we have 



v q 



e™d§dl (e- z ~ z f) = (-ir + y.q\ £ £ 



-i) J+fc H p _ M _ k (z,z) Q0k M 



3=0 fc= 



o J-k ] - {p-j)Kq-k)\ 



(2.11) 
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Proof. By applying repetitively the Leibnitz formula, it follows 



P Q 



j=0 k=0 



But from (2.2), i.e., e zz d^ (z n e 



[—l) m H mjn (z,z), we deduce 



v i 



^(^7) = EE(- 1VI 

j=0 fc=0 



This completes the proof. 



□ 



Remark 2.4. For f(z) = e~^ Z//2 with ( 6 C, £/ie right hand side of (2.11) leads to the 
generalized complex polynomials of Hermite type studied in [TT] and suggested by a special 
magnetic Schrodinger operator. 



Remark 2.5. Using (2.7) together with the established fact (—d z + zf {e zz g^j - 
we get 

{-ds + zf (-d z + zf (f) = (-l)^e 2 ^ (e-*V) • 
Therefore, the operational formula (2.11) can be reworded as 

<-* + " f + ^ (/) = ^ g £ ,1 b-^J^ cfl. 

Moreover, it is interesting to observe that H p ^ q (z,z) can be realized also as 

H m (z, z) = {-ds + zf (-d z + zf ■ (1). 



■lfe zz (of 



(2.12) 



(2.13) 



In the next section, we are concerned with some properties, for the complex Hermite 
polynomials, that follow directly from the previous operational representations, including 
( pE5| ), gjj and its variant Q as well as fl2~l~2] ). 

3. Related properties 

3.1. Recurrence formulae. Note for instance that, since d z and (— d z + zf commute, 
we get 

d z H M (z,z) = d z {-d z + zf (z p ) = {-d z + zfd z (z p ) = pH p ^, q (z,z). (3.1) 
Similarly, we have 

dgH M (z,z) = d z {-d z + zf{z«) = {-d z + zfd z {z^) = qH p>q ^(z,z). (3.2) 
Thus, we see 

(-d z + z) c\H m = qH m , (3.3) 

that is H pq are eigenf unctions of the second order differential operator of Laplacian type 
A := (— d z + z) d z . Furthermore, the polynomials H p q obey the recursion relations 

a) H„ + -i_„ = zH„ n — c\H„,„. (3.4) 

(3.5) 

l p,g+i — zHp^ — d z H PA . (3.6) 



lp+l,q ^ ±J -p,q lJ z- L - L p,q- 

b) H p+ i q = zH pq — qH pq _i. 
a) H n 



b') H, 



p,q+l 



zH PA - pH, t 



p-l,q- 



(3.7) 



Indeed, using (2.5), it follows 

H p+1 , q (z,z) = (-d^+z) 1+p (z«) 



-d z + z) H Ptq (z, z) = -chH m + zH m . 



Substitution of (3.2) in the previous equality gives rise to (3.5). Similarly (2.4) yields 
(3.6) and so (3.7). Recurrence formula (3.6) (resp. (3.7)) can also be obtained from 
(3.4) (resp. (3.5)) by conjugation since the polynomials H M (z,z) satisfy the following 
symmetry relation with respect to index permutation H Ptq (z,z) = H qp (z,z). 



3.2. Quadratic recurrence formulae. By writing (2.8) for the special case of 
/ = H p s , we get the following first variant of Nielsen's identity 

Hp,q->rS\Zi z , 



p ' q ' h kl (v-w (?-*)• ' 



(3.* 



This follows using H P!q+s 



-d z + z) q (H PtS ) combined with the fact that d k (H PjS 



pi 



H p _k jS for k < p. Furthermore, by considering the particular case / = z p+s in 



( 2.9[ ), we check that 



(-d z + zy(z? +s ) = q\J2 



k=0 



-l) k H M _ k {z,z) k 
k\ (q-k)\ A 



Whence, the second variant of the Nielsen identity reads 



H, 



p+s,q\Z } Zj 



min(q,s) f l^k z s~k jj 



, , IV ^ - L± p,q— k \Zj Z) 

q - S - to k\ (s-k)\ (q-k)\ 



(3.9) 



A third variant of Nielsen's identity can be established. More precisely, we have the 
following proposition. 



Proposition 3.1. We have 



H. 



i — ■ 

H p+S>q (z,z) = V2' 



■p+s-q 



s.k 



k=0 



z z 



H, 



p,q—k 



Z Z 



fc! 



(q-k)\ 



(3.10) 



Proof. The proof of (3.10) can be handled by writing e zz (— d z + z) q (z p+s e zz ) in two 
manners. Indeed, the first one follows by taking / = z p+s e~ zz in (2.7). In fact, 



e zz (-d z + zf (z p+s e~ zz ) = (-lfe 2zz d q z (z p+s e~ 2zz 

= s/2 q ~ P ~ S H p+s>q (V2z, V2z) 



upon an appropriate change of variable. The second one follows from (]2.9|), to wit 



e zz (-d z + zf (z p+s e- zz ) = q\e zz £ 



■l) k H m _k{z, z) a kr„s-zz\ 



fc=0 



, =0 k\ (q-k)\ 

"■ H s<k (z,z) H Ptq _ k (z,z 



d K Jz s e~ 



{q-k)\ 
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Thence, we have proved that 



H p+s>g (V2z, V2z) = \\ £ 



H s ,k(z, z ) H p q _ k (z, z 



k=0 



k\ 



(q-k)\ 



which gives rise to (3.10) by replacing z by z/y/2. □ 

We have called here ( 3.8[ ) (resp. ( |3.9 )) by the first (resp. second) variant of Nielsen 
identity, since we can write it as a weighted sum of a product of the same polynomials, 
according to the fact that z q ~ k = H Qq _ k (z,z) (resp. z s ~ k = H s _ k0 (z,z)). We reserve the 
appellation Nielsen's identity to the following 

Proposition 3.2. We have the following quadratic recurrence formula 

min(p,n) min(q,m) (_-,\j+k fx / x rr ( v v\ 

H p+m , q+n[ z,z)-p.q.m.n. ^ ^ {p _ j)l{q _ k)l {m _ k)l{n _ j)V 

(3.11) 



Proof. By taking / = H m ^ n (z,z) in (2.12), we obtain 

p « i-iy+ k H p ^ k ( z ,z) 



-d- z + zf (-d z + zf (H m , n (z, z)) = p\q\ E E ' -L " P ~^T'ZA d k AH m ,n{z, z)). 



Therefore the result (3.11 ) follows upon making use of 

+ zf {-d z + z) q {H m>n (z, z)) = H p+miq+n (z, z) 
together with the fact that 

4d k z (H m , n (z,z)) = rn\n\ Hm - k ^% 

[m — k)'\n — j)\ 

for k < m, j < n and d^d k (H mn (z,z)) = otherwise. 



□ 



As a special case of (3.11) when taking m = q and n = p (keeping in mind the fact 
that H mtn (z,z) = H n)Tn (z,z)), we state the following 

Corollary 3.3. We have 

-\y +k \H p ^ k ( z ,z)\ 2 



p q 



= £ m ((p _ m _ kyr 

3.3. Generating functions. 

Proposition 3.4. We have the following generating functions: 

+ °° yP 

«) Y,- H P,q(^) = (z-u) q e UZ - 

p=0 P- 

+°° v q 

<0 Y,^H m (z,z) = (z-vye™. 

q=0 I' 



(3.12) 



b) E 

p,q=> 

These equalities are valid for all u, v, z complex. 



+°° u p v q 

__ff ( 7 -~\ — p uz+vz-uv 

Pt q=o v- q- 



(3.13) 
(3.14) 
(3.15) 



Proof. The generating function a) is in fact the conjugate counterpart of a'). To prove 
a'), we begin by writing 



^0 ? ! 



-vd z +vz ^j?^ _ & vz e —vdz^j\ 

Now, since e^te*) = £ (f )(-w) fc ^" fc , it follows 



+00 v g 



(_<;) = (z - W )V*. 



This is exactly (3.14) from which we deduce 



HOC u p v q 



p,q=0 



00 ii p 

'•■ \ (z — v) p = e uz+vz ~ uv 



E 

p,q=0 P- 



□ 



Corollary 3.5. We have 



E '^M"- H P' 



^,2: J = e 



Furthermore, for every fixed integers m,n> 0, we have the identities 



-°°z p 



00 z q 



E -jH PtTl (z, z) = and ^ —H m:q (z, z) = 0. 



9=0 V 



According to the referee suggestion, it would be more interesting, in obtaining repro- 
ducing kernels using coherent states, if one can provide the sum 



+ ' x H Ptq (z,z)H P!q (w,w) 



E 

9=0 



p\q\ 



The following result gives a closed form of such sum. 

Proposition 3.6. For every z, w G C and fixed positive integers p, m, we have 

(-1) 



g H p>q (z, z)H m>q (w, w) 

9=0 



m\q\ 



Hp <rn (z -w,z- w)e % 



ml 



(3.16) 



Proof. Multiplication of both sides of (3.14) by e ww and next different at ing the obtained 
formula m-times with respect to w infers 



+00 1 

Y,-H p , q (z,z)d?(w q e 



ww \ 3m 



g =o y- 



d? (z - w) p e 



P e w{z-w) 



Now, in view of (2.3), i.e., <9™ (w q e 



'-l) m H qtm (w,w)e- wiD , we obtain 



+00 1 

E -^-M M { Z) z)H q ^{w,w) = ^L- e »»d? ((z - wfe^A . (3.17) 
^0 q\m\ ml v 7 
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The right hand side of the previous equation can be rewritten as 



-1)' 



ml 



e ww d m n z _ w y e w(z-w) 



■iy 



ml 

-l) m 

— j— « 
ml 

I2.:il (-1)' 



e (z-w)(z-w) Qin h^ z _ w y e -(z-w)(z-w) 

(-l) m e^<9™ (£ P e~^] 



i^=z—w 



ml 



-H p , m (z -w,z- w)e l 



Replacing this in (3.17) leads to the desired result. 



□ 



Remark 3.7. As a particular case of (3.16), when taking p = m, we get 



g H M (z,z)H (w,w) = _ 



(3.18) 



since in this case H PtP (£,£) = (— l) p p\iFi(— p; 1; |£| 2 ), see [H], where iFx(a;b;x) stands 
for the confluent hypergeometric function. Furthermore, by rewriting the sum in the left 

+oo p—l +oo 



hand side of (3.18) as = Yl + Yl > an d next taking z = w and n = q — p for q > p, 

q=0 q=0 q=p 

one deduces the following 



>p \H p+n ,p(z, z)\ 2 _ ( — 1 ) 

(P + n)\p\ p\ 



p _ p-i 
e z ~ z 



^2 \ h p,i( z ' z )\ 2 



q=0 



plql 



(3.19) 



Proposition 3.8. We have the following summation formula: 



+oo 



Hp,q{z, z)H PtS (u, u)Hq tr (v, v) _ ^(u+l) e ~(u+v+uv+l) 



p,q,r,s=0 



piqirlsl 



(3.20) 



Proof. We make use of (2.5) to write the left hand side of (3.20) as 



+oo 



p,q,r,s=0 y-'i-' - b - 



+oo i 

r!s! 



^° u p v q 



r,s=0 



\p,q=0 



pi ql 



By applying (3.15), we obtain 



+00 y 

p,q,r,s=0 P-'i-' 

+00 



uz+vz—uv 



+00 1 

E T7 (-^ + ^) r (-^ + w) s (e 

r,s=0 r - S - 

T \(~d v + VY ( ^ - e uz+v-z-uv 
+°° 1 



r=0 ' ■ 



r=0 



■1 + z-u) 

: e~ ~e~ 

v ! 



z(u+l) z(v+l) ^—(u+v+uv+l) 



□ 



Remark 3.9. T/ie special case of r = s = (m (3.20 ) ^ leads to the generating relation 



fl3.15[ ). 

3.4. Runge's addition formula. 

Proposition 3.10. We have the following addition formula 

H» (z + W , ~z + u) = pi* fi) (P+ " /2 E E ^ ^-* ( ^ ^ 



27 ^ 0fc = (q-k)\(p-j) ] - 

(3.21) 

Pnoo/. Set A 2 := -^=y + (^2z) and A z := + (V2z), and note that 

AfA n z ■ (1) = H m , n (V2z, V2z). (3.22) 
Note also that H p q (z + w,z + w) can be written in terms of A z and A w as 

H Ptq (z + w,z + w)= (--J— + (z + w)) (- d +(J+w)\ .(1) 

\ a{z + w) j \ d[z + w) J 



(2) {(A, + A U ) P (^ + ^) , }.(1). 



Now, since the involved differential operators A z , A^-, A w and A^ commute, the binomial 
formula yields 



h m (z + ^ = q (p+9)/2 f££ Q Q .(i) 

-(i)^f£e)(8(^)(^-« 



Thence, we obtain the asserted result according to (3.22). □ 
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Remark 3.11. An alternative way to prove (3.21 ) is to write H p q (z + w , z + w) in terms 
of the operator A z as 

l\(fH-?)/2 P 



H p , q (z + w,z + w) = (-) vr (4 + A w y .(( V2zy(V2wf->) 



and next applying the binomial formula, keeping in mind that the operators A z and A w 
commute and satisfy 

A r z (z s ) = 2~ s / 2 H s , r (V2z, y/2z). 
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